J /ip production in an equilibrating partonic system 



Xiao-Ming Xu 

Institute for Nuclear Theory, University of Washington, Box 351550, 

Seattle, WA 98195 
and 

Nuclear Physics Division, Shanghai Institute of Nuclear Research 
Chinese Academy of Sciences, P.O.Box 800204, Shanghai 201800, China 

Abstract 

Any color singlet or octet cc pair is created at short distances and then expands 
to a full size of J/tp. Such a dynamical evolution process is included here in cal- 
culations for the J/ip number distribution as a function of transverse momentum 
and rapidity in central Au-Au collisions at both RHIC and LHC energies. The cc 
pairs are produced in the initial collision and in the partonic system during the 
prethermal and thermal stages through the partonic channels ab — > cc^ s+l Lj] and 
ab — > cc[ 2S+1 Lj]x, and then they dissociate in the latter two stages. Dissociation of 
cc in the medium occurs via two reactions: (a) color singlet cc plus a gluon turns to 
color octet cc , (b) color octet cc plus a gluon persists as color octet. There are mod- 
est yields of cc in the prethermal stage at RHIC energy and through the reactions 
ab — > cc[ 2S+1 Lj] at LHC energy for partons with large average momentum in the 
prethermal stage at both collider energies and in the thermal stage at LHC energy. 
Production from the partonic system competes with the suppression of the initial 
yield in the deconfined medium. Consequently, a bulge within —1.5 < y < 1.5 has 
been found for the J ftp number distribution and the ratio of J /ip number distribu- 
tions for Au-Au collisions to nucleon-nucleon collisions. This bulge is caused by the 
partonic system and is thus an indicator of a deconfined partonic medium. Based 
on this result we suggest the rapidity region worth measuring in future experiments 
at RHIC and LHC to be -3 < y < 3. 

PACS codes: 24.85. +p, 12.38.Mh, 25.75.Dw, 25.75. Gz 

Keywords: Ultrarelativistic nucleus-nucleus collisions, Equilibrating partonic 
system, J/ip number distribution, Survival probability 
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1. Introduction 

A hot deconfined medium favors the dissociation of J/ip since enough hard gluons 
can overcome the large energy gap between the J/ip and a continuum state of cc [1]. 
Models based on perturbative QCD have shown that a dense partonic system can be pro- 
duced in central Au-Au collisions at RHIC and LHC energies [2-6] and then evolve toward 
thermal equilibrium and likely chemical equilibrium [7-11]. Such parton plasmas will be 
searched for soon in experiments at Brookhaven National Laboratory Relativistic Heavy 
Ion Collider (RHIC). The J/ip suppression has been taken as a thermometer to identify 
the evolution history of a parton plasma by showing transverse momentum dependence 
of the survival probability in the central rapidity region [12]. 

Charmonium melting inside a hot medium, which leads to J/ip suppression, was 
proposed by Matsui and Satz to probe the existence of the quark-gluon plasma [13]. 
Before the complete formation of charmonium is achieved, a pre-resonant cc is expanding 
from a collision point. Dominance of the color octet plus a collinear gluon configuration 
in the pre-resonance state [14] may account for the same suppression of ip' and J/ip 
production in proton-nucleus collisions [15]. The growth of the color octet configuration 
and its interaction with nucleons along its trajectory in a nucleus are essential ingredients 
in explaining measured J / ip production cross sections. In addition, the importance of color 
octet configurations has been verified in pp collisions at center-of-mass energy y/s =1.8 
TeV with the CDF detector at Fermilab [16]. Theoretically, the color-octet production 
at short distances and its evolution into physical resonances has been well formulated in 
nonrelativistic QCD [17]. At the collision energies of RHIC and LHC we can reasonably 
expect considerable contributions from the color octet mechanism. 

The evolution of ultrarelativistic nucleus-nucleus collisions, e.g. central Au-Au colli- 
sions at both RHIC and LHC energies, has been divided into three stages in Refs. [9,18,19]: 
(a) an initial collision where a parton gas is produced; (b) a prethermal stage where elas- 
tic scatterings among partons lead to local momentum isotropy [20]; (c) a thermal stage 
where parton numbers increase until freeze-out. The term 'partonic system' refers to 
the assembly of partons in the prethermal and thermal stages. The parton plasma only 
denotes the assembly of partons in the thermal stage. The cc pairs are produced in the 
initial collision, prethermal and thermal stages but disintegrate in the latter two stages. In 
order to understand and make predictions for J/ip yields of RHIC and LHC experiments, 
the following physical processes are taken into account, (a) In the initial collision, cc pairs 
are produced in hard and semihard scatterings between partons from incoming nuclei by 
2-^2 processes which start at order a 3 s through the partonic channels ab — > cc[ 2S+1 Lj]x. 
In the prethermal and thermal stages, cc pairs can also be produced in 2 — > 1 collisions 
which start at order a 2 via the partonic channels ab — + cc[ 2S+1 Lj] since partons in the 
deconfined medium have large transverse momenta, (b) The cc produced at short distance 
is in a color singlet (cc)i or color octet (cc) 8 configuration which has a certain probability 
to evolve nonperturbatively into a color singlet state. This J/ip production process is 
formulated in nonrelativistic QCD. (c) Since the color octet to singlet transition of (cc)s 
takes time, gluons in the partonic medium couple to the color octet state and destroy this 
transition process. Normally, dissociation cross sections for g + cc — > (cc) 8 depend on the 
pair size. Expansion of the cc from a collision point to a full J/ip size has to be taken into 
account. 
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A physical resonance formed by a cc pair may be one of J/ip, Xcj, ip' and others. 
Since the radiative transition from a higher charmonium state to the J/ip takes a much 
longer time, the transition of such a state with nonzero pt takes place outside the partonic 
system. Since the Fermilab Tevatron experiments have been able to separate direct J/^'s 
from those produced in radiative Xcj decays [16], in this work we assume that the direct 
J ftp production can also be extracted in heavy ion measurements. If Xcj and tp' are 
considered, suppression factors for Xcj and ip' in a deconfmed medium are included in 
prompt J/ip production. The identification of J/ip suppression in the medium becomes 
impossible for any prompt J/ip production data. Therefore, no contributions from higher 
charmonium states are taken into account in this work. 

The purpose of this work is to study the dependence of the J/ip survival probability 
and number distributions produced in central Au-Au collisions at RHIC and LHC energies 
on the transverse momentum and also rapidity which will be measured in RHIC experi- 
ments [21]. The J /ip number distributions corresponding to production of cc in the initial 
collision are given in Section 2. Since nuclear shadowing has been shown to influence J/ip 
production in proton-nucleus collisions [22], the nuclear modification of parton distribu- 
tions is considered. The J/ip number distributions due to cc production in the prethermal 
and thermal stages are given in Sections 3 and 4. Section 5 contains dissociation cross 
sections for gluon-(cc)i and gluon-(cc)s. Numerical results for nucleon-cc cross sections, 
J /if) number distributions and four ratios including survival probability are presented in 
Section 6. Conclusions are summarized in the final section. 



2. Initial production of cc 



Intrinsic transverse momenta of partons inside a nucleon result in the production of 
J/ip with typical momenta comparable to the QCD scale via 2 — > 1 partonic scattering 
processes [23]. Since we want to study J/ip productions with p? > 2 GeV, contributions 
from 2 — > 1 partonic reactions are not considered in the initial nucleon-nucleon colli- 
sion. The effect of intrinsic transverse momentum smearing is rather modest for large 
transverse momentum J/ip data from the Tevatron [24]. Upon omission of the intrinsic 
transverse momentum, differential cross section for J/ip production in nucleon-nucleon 
collision resulting only from 2 — » 2 partonic processes is given as 

d 3 <7 



dydy x dp ± 



2RL J2 X aXbfa/N(x a )f b/N (x b ) 
abx 

(i) dt 

+E§(^-^1 25+1 4 8) ]x-W)] (i) 

(8) 0,1 

where the summation J2 a bx is over partons labeled by a,b, x, f° r all possible color- 
singlet states and J2 (8) f° r all possible color-octet states. Here ^ denotes the partonic 
differential cross section for producing a cc[ 2S+1 Lj] and evolving to a J/^ with spectro- 
scopic notation for quantum numbers and superscripts for singlet and octet [23, 25], and 
fa/N is the parton distribution function of the species a in a free nucleon. The longitudinal 
momentum fractions carried by initial partons, x a and x b , are related to rapidities of cc 
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and x, y and y x , by 

x a = -^=(m±e y + p±e y *), x b = -^-=(m±e~ y + p±e~ yx ) 

V s V s 

where y/s, p± and are the center-of-mass energy of nucleon-nucleon collision, trans- 
verse momentum and transverse mass of the J/ip. The conditions x a < 1 and Xb < 1 
restrict y x to a region of 

, y/s — m±e~ y , y/s — m±e y 
- In < y x < In 

P± P± 

These 2 — > 2 processes at order a 3 , gg — > cc[ 2S,+1 Lj]g, gg — > cc[ 2S,+1 Lj]g, gg — > cc[ 2S,+1 Lj]g 
and gg — > cc[ 25+1 Lj]g, start in initial nucleus-nucleus collisions and proceed with the 
expansion of the heavy pair. While a cc propagates inside a prethermal or thermal partonic 
system, gluons hit and excite it to continuum states. Let a _ Il3c (i), be the cross section for 

gcc[l ibj ] 

g + (cc)[l 3 S[ 1] ] -> (cc) 8 , (7^5(8)] for g + (cc)[S {8) ] -> (cc) 8 and a^s), for (cc)[P (8) ] -> 
(cc) 8 respectively. The cross sections are calculated in Section 5. The probability for 
dissociation of a small-size cc into a free state relies on the relative velocity between the 
gluon and cc, v re i, and gluon number densities in the prethermal and thermal stages, n g (x) 
and n g (r), respectively. Here the variables x and r are individually space-time coordinates 
and proper time. In the prethermal stage, parton distributions depend on the correlation 
between momentum and space-time coordinates [18, 19]. The dependence of the gluon 
number density n g (x) on x characterizes the partonic system in nonequilibrium. In the 
thermal stage, thermal parton distributions can be approximated by Jiittner distributions 
where the temperature and parton fugacities depend only on the proper time [9, 18, 20]. As 
a consequence, the gluon number density is only a function of r. Including cc suppression 
in the partonic system, the finally-formed number distribution of J/ip resulting from cc 
pairs produced in the initial central A+B collision is given by 



mi 



dyd 2 pj 



/•In p± x fRA 

p± 

x a fa/A(x a , m 2 ± , r)x b fb/B(x b , m 2 , -f) 

/Tiso 
dr n g (x) < VreiVg^gm^k • u) > pre 0{d - V T At) 

dr'n g (r') < v rel a 3 sW {k ■ u) > the 6(d - V T At)] 

iso ^ 



+ ^(ab^cc[ 3 S ( * ) ]x^ jm 
dt 



r T iso 

exp[- / dr n g {x ) < v rel a gc£[S (s)](k ■ u) > pre 0(d - V T At) 

J TO 

- / dr'n g {r') < v ret a gcZ[s{&)] (k ■ u) > the 9(d - V T At)} 

J Tiso 



+ ^(ab^cc[ 1 si ) %^J/^) 
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f' T iso 

exp[- / dr n g (x ) < v re ia gc5[s(s)] (k ■ u) > pre 6(d - V T At) 

J TO 

f T f 

- / dT'n g (r') < v rel o g ^ [s ^)\{k ■ u) > the B(d - V T At)] 

J Ti er, 



rn ao 

exp[- / dr n g (x ) < ^ejCgccipC 8 )]^ • u) > pre 6(d - V T At) 

J TO 

f T f 

- / dT'n g (T') < Vreicrg^p^k ■ u) >the 9(d - V T At)]} 

(2) 

where f a /A is the parton distribution function of a nucleus, 

fa/ A (x, Q\ r) = T A (r)S a/A (x, f)f a/N {x, Q 2 ) (3) 

with the thickness function T\ and nuclear parton shadowing factor S a /A- Here, Ra 
is the nuclear radius. The symbols < ■ • • > pre and < • • • >the denote averages over 
gluon distributions in the prethermal and thermal stages, respectively. Along the track of 
nucleus-nucleus collisions, a deconfined partonic gas is produced from scatterings among 
primary partons at r , then reaches thermalization at r iso and finally freezes out at r/. 
Here, d is the shortest distance which a cc travels from a production point r to the surface 
of the partonic medium with transverse velocity Vp [12]. Suppose a cc is produced at a 
proper time r and a spatial rapidity rj. The time At for the partonic system to evolve to 
another proper time r' is 

(V\\ sinhr? — cosher + A /(sinhr7 — V\\ cosh i]) 2 T 2 + (1 — V, 2 )t' 2 

At = — V 2 " — (4) 

where Vm is the longitudinal component of the cc velocity. The disappearance of medium 
interactions on the cc is ensured by the step function 9 while this pair escapes from the 
partonic medium. 

3. Production of cc in the prethermal stage 

To order a 2 ., a cc in a color singlet state is produced only through gluon fusion 
gg — > cc[ 2S+1 L^]. For the ccpS*^], this fusion does not occur. In contrast, color octet 
states result from both channels gg — > cc[ 2S+1 L^] and qq — > cc[ 2S+1 L^]. Nevertheless, 
the number densities of quarks and antiquarks are so small that they are neglected in es- 
timating the production of cc in the prethermal stage where gluons dominate the partonic 
system. Four momenta of the two initial partons and final cc are denoted by k\ = (oui, hi), 
ki = ((V2,k2) andp= (E,p) = (m±coshy,p±,m±smhy). The differential production rate 
for gg — > cc[ 2S+1 L^] — > J/tjj in the prethermal stage is 

= f ^ fh fiWfa + k 2 - P )-g 2 G fg(k 1 ,X)fg(k 2 ,X) 

d s p 8(2vr) 5 i Wl u 2 V y)^GJg\ i, )Jg\ 2, > 

Y J \M{gg^cc[ 2S ^Lf]^J/^) 



(8) 



(5) 
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where go is the degeneracy factor for gluons and the f g (k, x) is the correlated phase-space 
distribution function given in Ref. [18]. The squared amplitudes | M. | 2 for cc in color 
singlet and color octet are calculated individually in Refs. [23, 25]. To order a 2 , the 
allowed color octet states are 1 Sq 8 ^ > and 3 Po,2 > through the gluon fusion channel. Taking 
into account the suppression of cc in the prethermal and thermal stages, the finally-formed 
number distribution of J/ip resulting from cc pairs produced through 2 — > 1 processes in 
the prethermal stage is given by 

dN 2 ^ 1 ttR 2 a r T is° k 2 1 

' rdrdr]d(f) kl dy kl ^-glf g (k 1 ,x)f g (k2,x) 



dyd 2 P± 16(2tt) 5 Jro ' ™ " M m 2 T 

{\M{gg^c^sf\^J^) \ 2 

dr'n g (x') < v re ia gc ^ s w](k ■ u) > pre 9(d - V T At) 

/Tf 
dr'ngir') < v re ia gc5[s{&)] (k ■ u) > the 9(d - V T At)] 
- iso 

+ \M{gg^cc[ 3 Pf ) ]^J/^) | 2 

dr'n g (x') < v re ia gc ^p(%)^{k ■ u) > pre 9(d - V T At) 

f T f 

- / dr'n g (r') < v re ia ^^(k ■ u) > the 9(d - V T At)]} 



(6) 



where (f) ki is the angle between k±i and pj_ for i — 1,2 and m c is the charm quark mass. 
The kinematic variables k±i, k± 2 , 4> k2 an d Uk 2 are expressed in terms of 



2m ' 



m± cosh(y - y kl ) - p ± cos(p kl 



k±2 = \lv±_ + k\ x - 2p ± k±i cos(p kl 

i k ±1 
sm0 fc2 = -- — sm0 fel 
k±2 

1 

sinh y k2 = - — (m± sinh y - k ±1 sinh y kl ) 
k±2 

To order a 3 , the differential production rate gets contributions from the processes 
gg -> cc[ 2S ^L^} g in the prethermal stage, 

d 3 A 2 ^ 2 1 r d 3 k 1 d 3 k 2 d 3 Px (4) 

E J = r-^r / —5 ( > (h + k 2 -P-Px) 

d A p 16(27t) s J uji U2 E x 



lg 2 G f 9 (ki,x)f 9 (k 2 ,x){Y: | M(gg - ccf^L^x -> J/rj,) 



2 (i) 
+ E I M(ffg - cc[ 2S+1 Lf]x -> J/^) | 2 } (7) 

(8) 

where p x = (E x ,p x ) = (p± x cosh y x ,p± x cos <p x ,p± x sin <p x ,p± x sinh y x ) is the four momen- 
tum of the massless parton x. Taking into account the suppression of cc in the prethermal 
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and thermal stages, the finally-formed number distribution of J/ ip resulting from cc pairs 
produced through 2^2 processes in the prethermal stage is given by 



dN 2 ^ 2 7iR 



— / TdTdr]p ±x dp ±x d(f) x dy x d(f)k 1 dy kl 

7T P J t () 



dyd 2 p ± 16(2tt) 8 J To 
2k 2 1 

—f- ^cfgih, x)f g (k 2 ,x) 
{\M{gg^cc[*S?]*^J/i>) \ 2 

dT'rig(x') < VreiPg^zsm^k ■ u) > pre 9(d - V T At) 

f T f 

- / drn g (T) < v re ia gc _, l3s (i)Ak ■ u) > the 9(d - V T At)] 

J T iso 1 

+ \M(gg^cc[ 3 s{%^J/^) \ 2 

dr n g (x ) < v rd a gc - c[s{s)] (k ■ u) > pre 9(d - V T At) 

dT'n g (T') < v rel a gc5[s(&)] (k ■ u) > the 9(d - V T At)} 

iso 

+ \M(gg^cc[ 1 S^]x^J/^) | 2 

dr n g {x ) < v re ia gc5[S (8)](k ■ u) > pre 9{d - V T At) 

- / dT'n g (r') < v rd a g C c[ S (»)](k ■ u) > the 9(d - V T At)} 

>>Tiso 

+ \M{gg^c<fpf ] ^^J^) | 2 

/7~iso 
dr'rigix') < Vreia^yp^ik ■ u) > pre 9(d - V T At) 

dr'ngir') < v rel a gc5[pm] (k ■ u) > the 9(d - V T At)]} 

iso 

where s — (ki + k 2 ) 2 and some kinematic variables are given by 

k±i = {Am 2 + 2m_Lj9_L x cosh(t/ - y x ) - 2p ± pj_ x cos0 x }/ 
{2[m_LCOsh(y - y kl ) +p± x cosh(y x - y kl ) -p±coscj) kl -Rl x cos(0 x - (p kl )]} 

k 2 ±2 = m\ + p\ x + 2m ± p ±x cosh(y - y x ) + k\ 1 

-2A;_ L i[m_ L cosh(y - y kl ) + p Lx cosh(y x - y kl )] 

sinh y k2 = j^—[m ± sinh y + p ±x sinh y x - k±i sinh y kl ] 

The J/tp number distribution resulting from cc pairs produced in the prethermal 
stage becomes 

rlN dN 2 ^ 1 dN 2 ~* 2 

ii±v pre "'-''pre _|_ " pre /g\ 

dyd 2 p± dyd 2 p± dyd 2 p± 
4. Production of cc in the thermal stage 



(8) 



In the thermal stage, parton distributions are approximated by thermal phase-space 
distributions fi(k; T, Aj) in which the temperature T and nonequilibrium fugacities \ are 
functions of the proper time r [9, 18]. While the partonic system evolves, quark and 
antiquark number densities increase. To order a 2 , both gg — > cc[ 2S+1 L^] — > J/ip and 
qq — > cc[ 2S+1 L^] —> J/ip contribute to the J/ip number distribution in the thermal stage 

dN 2 ^ 1 tiR 2 a m k 2 ±l 

{\ay t (hlT,X,)f t (h;T,X,) I M(gg -. ccfsSP} - J/« | 2 

dr'n g (r') < v re i<T g< ^ S {a)](k ■ u) > the 6(d - V T At)\ 
+ l -g 2 G f g {k 1 - T, X g )f g (k 2 ; T, A ff ) | M(gg -> ccfpf] - J/^) | 2 
exp[- y / dr'n g (T') < v re io gC c[P(8)](k ■ u) > the 6(d - V T At)] 
+g q9q f q (k 1 ; T, X q )f q (k 2 ; T, A ? -) | M(qq - cc[ 3 ^ 8) ] - J/V) | 2 
dr'n g {r') < v re i<T g c%sW](k ' u ) >the d{d - V T At)}} 

(10) 

where g q and g q are the degeneracy factors for quarks and antiquarks, respectively. In the 
channel of quark-antiquark annihilation, only the squared amplitude for 3 S[ 8 ^ does not 
vanish. 

All lowest-order 2-^2 reactions gg — > cc[ 2S+1 Lj]g, qq — > cc[ 2,s+1 Lj]g, gq — > 
cc[ 2S+1 Lj]g and — > cc[ 25+1 L j]g contribute to the J ftp number distribution in the ther- 
mal stage 

dN 2 ^ 2 ttR 2 a n-f ,,,,,, 2 ^ii 

l - ^ — = iwo \s / rdr dVP± x dp± x d(p x dy x d(j) kl dy kl — 
dyd 2 p± 16(27r) 8 J nso s 



\fa{ki\ T, \ a )f b (k 2 ] T, \ b )g a g b 

ab 

{J2 \M(ab^ cc[ 25+1 L? } ]x -> J/V) 
(i) 



dr'n g (r') < v rd a g ^ s+lL {.i)}{k ■ u) > the 9(d - V T At)} 
+ I M(ab -> cc[ 2S+1 Lf]x -> J/^) | 2 

(8) 

dT f n g (r') < w re /o- 9cE[i ( 8 )](A; • «) > tfte 0(d - V-rAt)]} 

(11) 

The J/ ip number distribution resulting from cc pairs produced in the thermal stage 
becomes 

dN the _ dN 2 ^ 1 | dNfc* 



dyd 2 p_i_ dyd 2 p± dyd 2 p 



_L 
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5. Gluon-cc dissociation cross sections 



A dissociation cross section of a full-size J/ip induced by a gluon is given in Refs. [1, 
26]. Since an initially-created cc has a radius of about r = and proceeds by expanding 
to a full-size object, the dissociation cross section of cc by a gluon has a size dependence. 
By this we mean the dissociation of cc into free states via this process g + cc — > (cc)g. 
Cross sections are calculated with chromoelectric dipole coupling between gluon and cc in 
the procedure for gluon- J/vf) dissociation in Ref. [26]. The wave function of an expanding 
cc is needed for this purpose, but it has not been investigated in the partonic medium even 
though some attempts have been made in studies of the color transparency phenomenon 
[27]. We proceed with the construction of wave functions in a simple one-gluon-exchange 
potential model. 

In a parton plasma, the internal motion of J/ip is obtained [12] from the attractive 
Coulomb potential, V\ = — g s 2 /37rr. The quantum-mechanical interpretation of the cc 
radius is V < r 2 >, the square root of the radius-square expectation value of the relative- 
motion wave function. For the IS color singlet, its wave function in momentum space 
normalized to the radius of cc[l 3 Sf^] is 

[*,](.V32^ ITT |2_ (13 ) 

where the variable a = \J < r 2 > /3 is the Bohr radius for a full-size J/ip. The velocity- 
square expectation value of the J/ip wave function is < v 2 >= 0.428. Then the radius 
of cc is assumed to expand according to V < r 2 > = V < v 2 >t + r . The gluon-cc[l 3 Si^] 
dissociation cross section is 

128g s 2 mfq 7 (Q° - e )^Q° 
a 9 cc[i^\ 9[m c a 2 (Q° _ €q) + i]6 ^ 

where Q° is the gluon energy, eo the binding energy of J/ip and g s the strong coupling 
constant. 

While the cc is in a color octet state, it is not a bound state but rather a scattering 
state. Its relative-motion wave function is determined by the repulsive potential V$ = 
g 2 /247rr. The radial part of the S wave function is 

S R (r) + iS T (r) = e^F(l + ir], 2, -2igr) (15) 

and the radial part of the P wave function is 

P R (r) + iPj(r) = qre iqr F(2 + irj, 4, -2iqr) (16) 

with q = m c V < v 2 > and rj = g 2 /247rV < v 2 >. The function F is the confluent hyper- 
geometric function. Wave functions in momentum space are obtained by performing a 
Fourier transform of the wave functions in space coordinates. Normalization constants of 
the momentum-space wave functions, C$ and Cp, are determined by fitting the cc radius. 
Dissociation cross sections of the S'-wave and P-wave color-octet states by a gluon are 

g 2 (m c Q ^- 5 - rb rb 



gWTp2 / / A A 3 3 

<W«] = 576. Cs L Jo 



ii(V^cQ ri)ji(Vm c g°r 2 )[^(r 1 )^(r 2 ) + SfrjSfo)] 



(17) 
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g s 2 (m e g°) L5 2 r 6 , rjVf 

^ (8) 1 = 576. Cp L I dndr2 ^ 

bo ( yj m c Q°ri ) j ( ^ m c Q°r 2 ) + 2 j 2 ( ^ m c Q n)j 2 ( ^ m c g°r 2 ) ] 

(18) 

where the j 5 Ji an d J2 are spherical Bessel functions. The b is determined so that the 
square root of the r 2 expectation value of the relative wave function in Eq. (15) or (16) is 
the color-octet radius. Relations b = 1.435V < r 2 > for S-wave and b = 1.3V < r 2 > for 
P-wave approximately hold for color-octet size less than normal hadron size. 



6. Numerical results and discussions 



Results for five aspects are presented in the following subsections. The first aspect 
is the nucleon-cc dissociation cross sections shown in the next subsection. The second one 
in Subsection 6.2 is J/ip number distributions versus transverse momentum at y = and 
rapidity at px — 4 GeV with nuclear effect on parton distributions and cc dissociation 
in the partonic system. The third one in Subsection 6.3 is to define and calculate four 
ratios including survival probability with y = or p T = 4 GeV at both RHIC and LHC 
energies. The fourth one is given in Subsection 6.4 to show J/ip number distributions 
without nuclear effect on parton distributions and cc dissociation in the partonic system. 
The fifth one concerns some uncertainties on the above results. 



6.1. Nucleon-cc dissociation cross sections 



In the parton model of the nucleon, the gluon is a dominant ingredient. Whereas 
the cross section for cc dissociated directly by a real gluon is of order a s , the cross section 
for the quark-cc dissociation through a virtual gluon is of order cc 2 . With the gluon-(cc)i 
cross section given in the last section, the nucleon-cc[l 3 S'j 1 ' ) ] cross section driven mainly 
by the gluon ingredient becomes 

a Ncc[l3s[V] = /(i) dx f S M X > Q 2 ) Vc[l3S«] (19) 

min 

where x^ in = ^ with being the proton momentum in the rest frame of the J/ip. The 
gluon distribution function f g / N is that Gliick-Reya-Vogt (GRV) result at leading order 
in Ref. [28]. The cross section is drawn in Fig. 1 to show the energy and renormalization- 
scale dependence while the cc radius is the J / ip radius in the attractive Coulomb poten- 
tial, rj/^ = 0.348 fm. In Fig. 1, gluon field operators are renormalized at three scales 
eo, V2eo,Q°, respectively. The coupling constant has the value a s = f^/^f- correspond- 
ing to the scale e while it varies for the other two scales. Values of the cross section 
at V^ = 10 GeV are a little lower than the nucleon- J ftp dissociation cross section ob- 
tained by the subtraction of quasi-elastic cross section in Ref. [29] from the total cross 
section given in Ref. [30]. In high-temperature hadronic matter or J ftp photoproduction 
reaction, a typical value of the center-of-mass energy for nucleon- J/ip (or preresonance) 
dissociation is around V^ = 6 GeV [1]. At this energy, Fig. 2 is drawn to show the size 
dependence of cr^^i)], with f g/N (x,Q 2 ) at Q 2 = e 2 ,. 
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For the S'-wave color octet the nucleon-ccfS 1 ®] cross section is 

°Ncc[S(8)] = / {8) dxf g/N (x, (g + Q°) 2 )Vc[5W] (2°) 

mm 

where a;£] n = ^ 3£R - For the 

P-wave color octet the nucleon-ccfP®] cross section is 

Vncc[pW] = / (8 ) dxf g/N (x, (q + g ) 2 ) o- flce[ p(8)] (21) 

Since the gluon momentum in a confining medium is bigger than the QCD scale [14] , the 
lowest value of x is set by Aq C d = 0.2322 GeV used in the leading order GRV parton 
distribution functions. Dependences of c Nc ^ S (8)^ and o"jv cg [p(8)] on the center-of-mass energy 
y/s are depicted in Fig. 3 while the size of (cc) 8 is the full size of J/ip. The dot-dashed 
line is obtained with the nucleon- J /if) cross section given by Eq. (24) in Ref. [1] where 
another gluon distribution function evaluated at Q 2 = is used. While the (cc) 8 has 
small momentum in a nucleus, the cross section for nucleon- (cc) 8 production is lower than 
the absorption cross section determined by Gerschel and Hiifner [31] or the two-gluon 
exchange result [32]. It was proposed by Kharzeev and Satz that the color octet plus a 
gluon configuration is a dominant component produced in the proton-nucleus collisions 
[14]. In fact, the present cross section for a nucleon and a bare (cc) 8 is one part of the 
nucleon-g(cc) 8 cross section. 

A (cc) 8 pair produced at a collision point expands before becoming color singlet to 
a size which may be larger or smaller than the full size of J/ip. We then show in Fig. 4 
the dependence of the nucleon-(cc) 8 cross section on the color-octet pair radius. 

We do not want to address proton-nucleus collisions in terms of nucleon-cc cross 
sections [33, 34] since only the gluon-cc cross sections are needed to study cc suppression 
in the prethermal and thermal stages. In proton-nucleus collisions, once a color-octet 
(cc) 8 pair is produced, it picks up a collinear gluon to form a colorless configuration [14]. 
However, in central Au-Au collisions at RHIC and LHC energies, the accompanying gluon 
scatters with other hard gluons in the dense partonic system and is driven away. Therefore 
the bare cc is the object that we want to study in the partonic system. 

6.2. J/ip number distributions with suppression 

J/ip number distributions versus transverse momentum at y = and rapidity at 
Pt = 4 GeV for central Au-Au collisions at RHIC energy ^Js = 200AGeV are calculated 
with respect to the initial collision, prethermal and thermal stages. Initial productions 
of cc are calculated with GRV parton distribution functions at renormalization scale \i = 
+ Am? c . Evolution of color-octet states 3 Si, 1 S and 3 Pj toward the J/tp is specified 

by nonperturbative matrix elements < O^^Si) >, < O^^So) > and < (9^( 3 P ) > 
in nonrelativistic QCD [17]. In the nonperturbative evolution, a gluon from the partonic 
system hits and prevents the color octet from color neutralizing via g + (cc)s — > (ce)$. 
This medium effect has been expressed by exponentials in Eqs. (2), (6), (8), (10) and 
(11). Therefore, the nonperturbative matrix elements are assumed to be invariant while 
the medium effect is factorized into exponential forms. Values of these matrix elements 
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are well determined by fitting the CDF measurements for pp collisions at y/s =1.8 TeV 
in Ref. [35], 

< Og^fSi) >= (1-12 ± 0.14) x l(T 2 GeV 3 (22) 

< O^CSo) > +^4 < O J /\ 3 P ) >= (3.90 ± 1.14) x 10- 2 GeV 3 (23) 

m c 

In pp collisions, differential cross sections of direct J/ip production depend on the com- 
bination of < Ol^^So) > and < C 8 J/V, ( 3 Po) >• However, since the dissociation cross 
section for the S'-wave color-octet state is different from that for the P-wave color-octet 
state, such a dependence on the combination is destroyed. In calculations, values are 
taken as follows, 

< O^CSq) >= 4 x 10- 2 GeV 3 , < O J 8 ^( 3 P ) >= x 10" 2 GeV 3 (24) 

35 

The value of < 0^( 3 Pq) > is positive at tree level and negative after renormalization 
[36]. Eq. (23) is still satisfied by the values in Eq. (24). These values of nonperturbative 
matrix elements are supposed to be universal for any center-of-mass energy y/s. 

Various contributions to the J/ip number distributions including initial collisions, 
prethermal and thermal stages, 2 — > 1 and 2^2 collisions, are drawn separately in Figs. 
5 and 6. The dashed curve resulting from cc production in the initial collision is obtained 
by calculating Eq. (2) where the nuclear parton shadowing factor is given in Ref. [4] 
throughout this subsection. The upper and lower dot-dashed curves resulting from cc 
production in the prethermal stage are obtained by individually calculating Eq. (6) for 
2 — > 1 collisions and Eq. (8) for 2 — * 2 collisions. The upper and lower dotted curves 
resulting from cc production in the thermal stage are obtained by calculating Eq. (10) 
for 2 — > 1 collisions and Eq. (11) for 2 — > 2 collisions, respectively. To exclude the effect 
of intrinsic transverse momentum smearing, only the region px > 2 GeV is considered. 
Consequently, no 2 — > 1 collisions contribute in the initial collision. The J/x/j number 
distribution resulting from the initial collision shown by the dashed line has a plateau 
similar to that in proton-proton collision [37]. Both Figs. 5 and 6 show that cc pairs 
produced from the thermal stage can be neglected compared to the initial production, but 
the contributions from the prethermal stage are important in the transverse momentum 
region 2GeV < px < 8GeV and rapidity region < y < 1.2. Productions of cc in the 
prethermal and thermal stages bulge up the J / ip number distribution shown by the solid 
line in this rapidity region. Nevertheless, the dot-dashed and dotted lines fall rapidly as 
the rapidity gets large. This bulging characterizes the formation of a deconfined medium 
because the medium has average momentum limited but big enough to produce extra cc 
and thus J ftp. The 2 — > 1 collisions in the partonic system have bigger contributions than 
the 2 — > 2 collisions. 

Each of Figs. 7 and 8 contains two sets of lines to show contributions from the 
color-singlet and color-octet pairs produced at short distance. Any set has a dashed line 
obtained from Eq. (2) for the initial collision, a dot-dashed line from Eq. (9) for the 
prethermal stage and a dotted line from Eq. (12) for the thermal stage, respectively. A 
line in the upper (lower) set for the color-octet (color- singlet) contributions stems from 
the terms for (cc)s ((cc)i) states. The color octet states dominate productions of J/ip at 
RHIC energy. However, the ratio of color-octet to color-singlet contributions shown by the 
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two solid lines at p T = 6 GeV is reduced from about 70 at CDF collider energy y/s =1.8 
TeV to about 40 at RHIC energy. Both contributions of color-singlet and color-octet 
states have similar dependence on transverse momentum and rapidity. 

Figs. 9 and 10 show transverse momentum and rapidity dependence of J/ip number 
distributions for central Au-Au collisions at LHC energy y/s = 5.5ATeV. A prominent 
feature is that the J ftp number produced from the thermal stage is comparable to that 
from the prethermal stage. Compared to the initial production, cc and J/ip produced 
through 2^2 reactions may be neglected. A bulge is observed on the plateau in the 
rapidity region < y < 1.5. Such a bulge can be taken as a signature for the existence 
of a parton plasma at the LHC energy. Figs. 11 and 12 depict contributions from the 
color singlet and color octet at LHC energy. The ratio of color-octet to color-singlet 
contributions shown by the two solid lines at pr = 6GeV reaches about 150. This indicates 
the color octet states become more important with the increase of \/s. 

6.3. Ratios including J/ip survival probability 

Nuclear shadowing results in a modification of gluon distribution functions inside 
a nucleus [38] and such a nuclear effect is represented by the shadowing factor S a /A in 
Eq. (3). If the S a /A = 1 for no shadowing, the dN 2 ~^ 2 /dyd 2 p± is proportional to the 
product of atomic masses of the two colliding nuclei. If S a /A ^ 1 and depends on the 
longitudinal momentum fraction x, the production of cc is reduced in the shadowing region 
and enhanced for the anti-shadowing region. Irrespective of interactions of J/ip with the 
partonic system, J/ip number distributions produced in the initial central A+B collision 
is obtained by putting all exponentials equal to 1 in Eq. (2), 

d%l {Sa/A) = Um A.- d n drr 



dyd 2 p_\ 



X afa/A(Xa, m\, r)x b f b/B (x b , Ui\, -f) 

{^(ab^cc[ 3 S{%^JM 
+ ^(ab^cc[ 3 S[ 8) ]x^ J/ij) 
+ ^ab^cc[ l Sj,%^ J/VO 
+ ^(a6-cc[ 3 pf]x^ J /i,)} 



(25) 



To characterize the influence of nuclear parton shadowing on the J/ip production from 
the initial collision, a ratio is defined as 

Here the S a /A ^ 1 from Ref. [4] applies throughout this subsection. 

The initially produced J / ip originates from the cc pairs produced in the initial colli- 
sion. Its dependence on the transverse momentum and rapidity is obtained by calculating 
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Eq. (25). Some cc pairs produced in the initial collision may dissociate by gluons from the 
partonic system. As a consequence, the J/if> number is reduced. The survival probability 
for the cc transiting into a J /if) is defined as the ratio 

dN 2 ^ 2 dN 2 ^ 2 

^ = _g_ ( ^^ 1)/ _^_ ( ^ /1) (27) 

We have calculated the J /if) number distributions produced in the prethermal and 
thermal stages in Subsection 6.2. The J /if) yield may be bigger than the reduced amount 
of initially produced J / if) due to the cc dissociation by gluons in the partonic system. The 
partonic system has two roles. One is to produce cc pairs and another is to dissociate cc 
pairs. To see the roles, a ratio is defined by 

R plasma = / dN i^j 2 / g ^ ]_) + dN Pre + dN the ^ / dN Q^ A ^ \) (28) 

dyd 2 p± a ^ A dyd 2 p±_ dyd 2 pj_ dyd 2 p± a ^ A 

To understand the nuclear effect on parton distributions and the roles of the partonic 
system, we need to compare J/ if) production in the central A+B collision with that in the 
nucleon-nucleon collision. To this end, a ratio is defined as 

p / dN irti 2 f g _L 1 \ | dN P re | dNthe ^ / dN 0~* 2 ( c I \ l r )Ok\ 

R = ( dy^K {Sa/A * 1} + dy^K + ^V V ^V ( ^ = 1} (29) 
which is also written as 

R = R im RP lasma (30) 

The ratios R mi , S plasma , RP lasma and R versus transverse momentum and rapidity 
are depicted as dashed, dotted, dot-dashed and solid lines, respectively, in Figs. 13 and 14 
for the RHIC energy and Figs. 15 and 16 for the LHC energy. In contrast to R mi < 1, the 
value of RP lasma i s larger than 1 for all transverse momenta in Fig. 13 and < y < 1.5 in 
Fig. 14 and 0.5 < y < 1.5 in Fig. 16. This results in prominent bulges on the solid lines 
of R in Figs. 14 and 16. In contrast, the survival probability S plasma shown by the dotted 
lines has no such bulge. Therefore, the bulges are present in Figs. 6 and 10 when the J /if) 
yield resulting from cc pairs produced in the partonic system overwhelms the reduced 
amount of initially produced J /if). We conclude that in the rapidity region < y < 1.5 
a bulge observed in the ratio R is an indicator for the existence of the partonic system. 
For R im < 1 and S plasma < 1, J /if) suppression arises from the nuclear parton shadowing 
found in HIJING and cc dissociation in the partonic system. 



6.4- J /if) number distributions with no suppression 



In Subsection 6.2, J /if) number distributions have been presented while the cc re- 
duction due to the nuclear parton shadowing in the initial collision and cc dissociation in 
the partonic system are taken into account. In this subsection, the suppression including 
both the reduction and dissociation is omitted in calculations of J / if) number distributions 
by setting to 1 all exponentials in Eqs. (2), (6), (8), (10) and (11). Figs. 17-20 depict 
these distributions versus transverse momentum at y = and rapidity at pt = 4 GeV at 
both RHIC and LHC energies. We are now ready to explain the dip within y < 1 in Fig. 
10. This dip disappears in Fig. 20 where suppression is not considered. Since R mi shown 
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by the dashed line in Fig. 16 is flat with respect to the rapidity y < 2 and S plasma shown 
by the dotted line has a steep rise in 0.5 < y < 1, the dip phenomenon is solely due to the 
cc dissociation in the partonic system. Such a dip phenomenon is not obvious but still 
can be observed in the prethermal and thermal stages when the J / ip number distributions 
with suppression are compared to those without suppression. The comparison is indicated 
in Fig. 21 for the prethermal stage and Fig. 22 for the thermal stage. The solid and 
dot-dashed lines for no suppression begin to fall from y = 0.5 to y — 1, but change to 
rising as shown by the dashed and dotted lines when the cc dissociation is switched on. 
This change occurs because of the steep rise of S plasma . A relatively weak dependence of 
gpiasma on ^ j g s ] lowri the dotted line in Fig. 15. The dip phenomenon thus cannot 
be observed in the pt dependence of J/tp number distributions. Referring back to Eqs. 
(2), (6), (8), (10) and (11), exponentials there have sensitive dependence on the rapidity 
in < y < 1.5. The dip is more obvious in the color singlet channel as shown by the 
lower solid, dashed, dot-dashed and dotted lines in Fig. 12. This is so because the cross 
section for the l^^-state dissociation has a narrower peak with respect to the incident 
gluon energy [12] than the color octet states. 

6.5. Uncertainties 

Since gluon shadowing in nuclei has not been studied experimentally, theoretical 
estimates of the nuclear gluon shadowing factor involve uncertainties. The nuclear parton 
shadowing factor found in HIJING [4] is a result of the assumption that there is no Q- 
dependence on the shadowing factor and the shadowing effect for gluons and quarks is the 
same. Nevertheless, the shadowing factor has been shown by Eskola et al. to evolve with 
momentum Q [39] . The difference between the latter and the former indicates uncertainty. 
The ratio R mi defined in Eq. (26) is calculated with Eskola et a/.'s parametrization [39] 
and results are depicted in Fig. 23 showing momentum dependence at y — and Fig. 
24 showing rapidity dependence at pt = 4 GeV. Compared to the dashed lines in Figs. 
13-16, the change of R tm at RHIC energy greater than 1 is prominent. This implies that 
the anti-shadowing effect of Eskola et a/.'s parametrization is quite important at RHIC 
energy. Measurements on R mi in RHIC experiments are needed to confirm this nuclear 
enhancement [21]. 

The ratio R mi is always flat within the rapidity region -1.5 < y < 1.5 for parametriza- 
tions given in HIJING and by Eskola et al, and the flatness seems to be independent of 
parametrizations. If the partonic system does not come into being, the ratio R is flat, too, 
since R = R im . If the partonic system dissociates cc pairs, the solid curve of R undergoes 
bulging, dipping and then bulging from y = —1.5 to y — 1.5. Any such twist of R in 
— 1.5 < y < 1.5 observed in experiments is nontrivial, because only a deconfined medium 
generates it. 

Upon inclusion of uncertainties on the formation and evolution of parton plasma 
arising from other factors, for instance, the dependence on the coupling constant a s 
[40] and transverse flow [41], the J/ip number distribution and the four ratios including 
survival probability will change. In the partonic system considered here gluons dominate 
the evolution and gluon-cc interactions break the pairs. In a system where quarks and 
antiquarks are abundant, interactions between quarks (antiquarks) and cc may account 
for a suppression of J/ip [42]. Additional suppression caused by energy loss of the initial 
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state has not been considered since there is a controversy on the influence of the energy- 
loss [43, 44]. Some uncertainties are expected to be fixed by upcoming experiments at 
RHIC. 

7. Conclusions 

We have studied J/ip production through both color-singlet and color-octet cc chan- 
nels with various stages of central Au-Au collisions at both RHIC and LHC energies. In 
addition to the scattering processes ab — > cc[ 2S+1 Lj]x, contributions of the reactions 
ab — > cc[ 2S+1 Lj] are also calculated in the prethermal stage and thermal stage. The effect 
of the medium on an expanding cc involves a gluon interacting with the cc to prevent it 
from a transition into a color singlet. Cross sections for g + cc — > (cc)g are calculated 
with internal wave functions of (cc)i in an attractive potential and (cc)s in a repulsive 
potential. Furthermore, nucleon-cc cross sections for color singlet, S- and P- wave color 
octets as a function of y/s or cc radius are evaluated by assuming that the nucleon domi- 
nantly contains gluons. Momentum and rapidity dependence of J/ ip number distribution 
with various contributions are calculated for central Au-Au collisions at both RHIC and 
LHC energies. Color octet contributions are one order of magnitude larger than the color 
singlet contributions. Yields of cc are large in the prethermal stage at RHIC energy and 
through the 2 — > 1 collisions ab — > cc[ 2S+1 Lj] at LHC energy. Since the partonic system 
offers fairly large amounts of cc, a bulge in < y < 1.5 at RHIC energy and 0.5 < y < 1.5 
at LHC energy can be observed in the rapidity dependence of the J/ip number distribu- 
tion and the ratio R of J/ip number distributions for Au-Au collisions to nucleon- nucleon 
collisions. Such a bulge is a signature for the existence of a deconfined partonic medium. 
We suggest that RHIC and LHC experiments measure J ftp number distributions and 
the ratio R in the rapidity region < y < 3 to observe a bulge. While the yield of cc 
from the medium is larger than the reduced amount of initial production in the medium, 
the ratio RP lasma is larger than 1. The competition between production and suppression 
determines the values of RP lasma ^ which relies on the evolution of parton number density 
and temperature of the partonic system [12]. A dip in the rapidity dependence of the J/ip 
number distributions at LHC energy may exist and this amounts to a suppression effect 
of cc in the partonic system. So far, we have obtained results and conclusions for positive 
rapidity. It is stressed that the same contents for negative rapidity can be obtained from 
the positive region by symmetry. 
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Figure 1: Solid, dashed and dot-dashed lines are nucleon-ccfl 3 ^ 1 ^ cross sections for 
fg/N(x,Q 2 ) evaluated at Q 2 = el,2el, (Q ) 2 , respectively. The (cc)i has the same size as 
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Figure 2: Cross section for nucleon-cc[l 3 5'} 1 ' ) ] at = 6 GeV as a function of the cc radius 
is calculated with f g /N(x,Q 2 ) evaluated at Q 2 = eg. 
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Figure 3: The solid and dashed lines are cross sections for nucleon-ccfS 1 ^] and nucleon- 
C c[p( 8 )] collisions as a function of y/s, respectively. The dot-dashed line is the nucleon-J/^ 
cross section calculated with Eq. (24) in Ref. [1]. The corresponding (cc)g and J/ip have 
the same radius. 
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Figure 4: The solid and dashed lines show radius dependence of cross sections for nucleon- 
(cc)s[S^] and nucleon-(cc)s[-P^ 8 - ) ] collisions, respectively. 
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Figure 5: J/ip number distributions versus transverse momentum at y = and RHIC 
energy with suppression. The dashed curve corresponds to cc production in the initial 
collision. The upper and lower dot-dashed (dotted) curves correspond to cc produced 
through 2 —>■ 1 and 2^2 reactions in the prethermal (thermal) stage, respectively. The 
solid curve is the sum of all contributions. 
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Figure 6: The same as Fig. 5, except for rapidity distribution at p? = 4 GeV. 
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Figure 7: J/ip number distributions versus transverse momentum at y = and RHIC 
energy with suppression. The upper and lower dashed (dot-dashed, dotted and solid) 
lines correspond to cc in color octet and color singlet, respectively, produced in the initial 
collision(prethermal stage, thermal stage and the all three stages). 
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Figure 9: J/ip number distributions versus transverse momentum at y = and LHC 
energy with suppression. The dashed curve corresponds to cc productions in the initial 
collision. The upper and lower dot-dashed (dotted) curves correspond to cc produced 
through 2 — >• 1 and 2^2 reactions in the prethermal (thermal) stage, respectively. The 
solid curve is the sum of all contributions. 
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Figure 10: The same as Fig. 9, except for rapidity distribution at pr = 4 GeV. 
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Figure 11: J/tp number distributions versus transverse momentum at y = and LHC 
energy with suppression. The upper and lower dashed (dot-dashed, dotted and solid) 
lines correspond to cc in color octet and color singlet, respectively, produced in the initial 
collision (prethermal stage, thermal stage and all three stages). 
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Figure 12: The same as Fig. 11, except for rapidity distribution at px = 4 GeV. 
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Figure 13: Ratios versus transverse momentum at y = and RHIC energy. The solid, 
dashed, dot-dashed and dotted lines are R, R in \ W lasma and S plasma , respectively. 
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Figure 14: The same as Fig. 13, except for rapidity distribution at = 4 GeV 
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Figure 15: Ratios versus transverse momentum at y = and LHC energy. The solid, 
dashed, dot-dashed and dotted lines are R, R in \ RP lasma and S plasma , respectively. 
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Figure 16: The same as Fig. 15, except for rapidity distribution at px = 4 GeV. 
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Figure 17: J/ip number distributions versus transverse momentum at y — and RHIC 
energy without suppression. The dashed curve corresponds to cc productions in the initial 
collision. The upper and lower dot-dashed (dotted) curves correspond to cc produced 
through 2 — > 1 and 2^2 reactions in the prethermal (thermal) stage, respectively. The 
solid curve is the sum of all contributions. 
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Figure 18: The same as Fig. 17, except for rapidity distribution at px = 4 GeV. 
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Figure 19: J/tp number distributions versus transverse momentum at y — and LHC 
energy without suppression. The dashed curve corresponds to cc productions in the initial 
collision. The upper and lower dot-dashed (dotted) curves correspond to cc produced 
through 2 — >• 1 and 2^2 reactions in the prethermal (thermal) stage, respectively. The 
solid curve is the sum of all contributions. 
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Figure 20: The same as Fig. 19, except for rapidity distribution at px = 4 GeV. 
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Figure 21: J/ip number distributions versus rapidity at = 4 GeV in the prethermal 
stage of LHC energy. The dashed and solid lines individually correspond to cc produced 
through 2 — ► 1 reactions with and without suppression. The dotted and dot-dashed lines 
through 2 — > 2 reactions with and without suppression, respectively. 




Figure 22: The same as Fig. 21, except for the thermal stage. 
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Figure 23: Ratio R mi versus transverse momentum at y = is calculated with Eskola's 
parametrization. 




Figure 24: The same as Fig. 23, except for rapidity distribution at px — 4 GeV. 
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